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INTRODUCTION 

A  recent  addition  to  the  tools  of  management  and  production  engi- 
neering is  the  technique  of  linear  programming.  The  purpose  of  this 
project,  a  result  of  a  survey  of  literature  in  this  field,  is  to  sum- 
marize the  fundamentals  of  linear  programming.  Toward  that  objective, 
the  following  phases  are  presented:  a  summary  of  the  mathematical  con- 
cepts involved,  three  methods  of  solution  and  their  characteristics, 
and  the  significance  of  the  dual  theorem.  The  general  form,  the  simplex 
method,  and  two  of  its  special  cases,  the  transportation  method  and  the 
assignment  method  are  presented.  Two  basic  references  for  the  mathe- 
matics of  linear  programming  are  references  (l)  and  (U). 

The  names  of  the  special  cases  should  not  be  taken  literally.  For 
example,  the  transportation  method  may  be  used  for  problems  with  no 
physical  transportation  involved.  Similarly,  the  assignment  method  may 
be  used  in  quality  control,  purchasing,  set  up  time,  and  such  other  di- 
verse problems.  In  the  development  of  the  special  cases,  the  type  of 
problem  was  essentially  the  same  as  the  name  implied,  but  in  recent 
years  the  methods  have  been  applied  to  a  variety  of  problems. 

The  general  form,  the  simplex  method,  was  developed  by  Dantzig." 
Most  of  the  work  in  the  field  of  linear  programming  has  been  conducted 
by  economists,  mathematicians,  and  consulting  firms.  However,  the  num- 
ber of  qualified  personnel  for  this  type  of  work  is  limited  compared 
to  the  needs  of  industry.  The  usefulness  of  mathematical  programming 
depends  to  a  large  extent  on  the  ingenuity  of  the  personnel  formulat- 
ing the  mathematical  model.  Many  management  and  technical  publications 
contain  reports  on  the  use  of  this  type  of  programming  but  are  neces- 
sarily brief  due  to  publication  requirements  and  other  factors. 
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It  should  be  noted  that  many  of  the  calculations  involved  in  lin- 
ear programming  can  be  handled  in  a  routine  fashion  utilizing  desk  cal- 
culators, but  a  large  scale  problem  is  better  solved  by  computers. 
Regarding  the  use  of  desk  calculators,  one  of  the  advantages  of  linear 
programming  is  that  the  mechanics  of  solution  are  simple  and  may  be 
handled  by  lesser  skilled  personnel  to  obtain  the  best  solution  to  a 
given  problem.  This  characteristic  should  be  borne  in  mind  as  an  aid 
in  freeing  engineering  and  management  personnel  for  more  demanding  work 
commensurate  with  their  training.  The  most  important  phase  of  the  pro- 
gramming solution  is  that  of  setting  up  the  initial  problem.  Once  this 
has  been  done,  the  mechanics  can  be  performed  by  clerical  personnel. 
Operations  research,  operations  analysis,  and  mathematical  pro- 
gramming have  received  wide  publicity  since  World  War  II.  The  majority 
of  the  research  and  application  was  conducted  during  World  War  II,  and 
management  is  now  in  the  process  of  utilizing  some  of  the  methods  in 
industry.  There  is  an  abundance  of  literature  concerning  the  above 
topics  and  no  attempt  will  be  made  to  describe  or  disciiss  those  fields. 
Essentially,  linear  programming  is  a  method  or  technique  within  the 
field  of  mathematical  programming,  and  attention  will  be  limited  to 
that  phase. 

Linear  programming  is  a  systematic  mathematical  approach  to  a 
given  problem  to  arrive  at  an  optimal  or  best  solution.  It  involves 
an  objective  and  a  set  of  restrictions  or  constraints.  The  procedure 
has  been  effectively  used  in  industrial  problems  where  a  large  number 
of  alternative  choices  were  available.  From  case  histories  published, 
it  may  be  said  that  the  usefulness  of  the  method  is  proportional  to 
the  size  and  complexity  of  the  situation.  Until  more  research  has  been 
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conducted,  some  industries  may  not  be  able  to  use  linear  programming  to 
its  full  capacity.  However,  by  examining  a  local  problem  for  applica- 
tion and  trying  to  formulate  the  problem  by  the  methods  of  linear  pro- 
gramming, possible  cost  and  time  savings  may  be  realized. 

One  of  the  limitations  of  the  method  is  the  requirement  of  linear- 
ity of  the  variables  involved.  Due  to  this  restriction,  all  variables 
such  as  time,  cost,  and  quantity  must  be  of  the  first  order.  That  is, 
the  ratio  of  input  to  output  must  be  a  fixed  ratio  and  independent  of 
the  level  of  production  or  process  involved.  For  example,  if  one  unit 
of  input  results  in  one  unit  of  output  then  ten  units  of  input  must  re- 
sult in  ten  units  of  output.  Dorfman^  has  summarized  the  basic  assump- 
tions of  linear  programming  as  follows: 

1.  The  productive  opportunities  of  an  economy  or  economic 
unit  are  defined  by  the  resources  and  the  productive 
processes  available  to  it.  The  quantities  of  at  least 
some  of  the  resources  are  finite  and  so  is  the  number 
of  productive  processes  available. 

2.  Any  productive  process  may  be  used  at  any  positive  level 
consistent  with  the  supply  of  resources  available.  The 
consumption  of  resources  and  the  output  of  products  is 
proportional  to  the  level  at  which  the  process  is  used. 

3«  Several  productive  processes  may  be  used  simultaneously 
if  the  supply  of  resources  is  adequate.  If  this  is  done, 
the  consumptions  of  the  individual  processes  used,  and 
the  output  of  products  is  the  sum  of  the  outputs  of  the 
individual  processes. 

Dorfman^  further  reduces  the  three  assumptions  to  the  following 

postulates: 

a.  Linearity.  By  definition,  in  linear  programming,  each 
process  is  characterized  by  certain  ratios  of  the  quan- 
tities of  the  inputs  to  each  other  and  to  the  quantities 
of  each  of  the  outputs.  These  ratios  are  defined  to  be 
constant  and  independent  of  the  extent  to  which  the  proc- 
ess is  used. 


b.  Divisibility.   It  is  assumed  that  any  process  can  be 
used  to  any  positive  extent  so  long  as  sufficient  re- 
sources are  available;  indivisibilities  and  "lumpiness" 
in  production  are  ignored. 

c.  Additivity.  It  is  assumed  that  two  or  more  processes 
can  be  used  simultaneously,  within  the  limitations  of 
available  resources,  and  that  if  this  is  done  the  quan- 
tities of  the  outputs  and  inputs  will  be  the  sums  of 
the  quantities  which  would  result  if  the  several  proc-  ( 
esses  were  used  individually. 

d.  Finiteness.  It  is  assumed  that  the  number  of  processes 
available  is  finite. 

The  terms:  resources,  products,  and  production  processes  have  very 
broad  meanings  as  used  above  and  are  not  limited  to  the  concepts  often 
used  in  industry. 

Programming  in  this  sense  may  be  considered  as  arriving  at  the  solu- 
tion systematically;  the  term  stems  from  the  "activity  analysis"  approach 
of  reference  (ii).  The  objective  of  a  problem  may  be  to  minimize  or  maxi- 
mize some  function  such  as  time,  cost,  or  material.  The  constraints  may 
be  set  by  machine  time  available,  raw  materials,  capital  to  be  invested, 
time  in  transportation,  quantity  to  be  produced,  or  other  restrictions. 

Initially,  the  particular  problem  is  formulated  mathematically  into  a 
set  of  inequalities  with  more  variables  present  than  the  number  of  equations. 
If  it  is  possible  to  set  up  as  many  equations  as  variables  present,  then  the 
obvious  method  is  a  solution  by  simultaneous  equations.  However,  in  such  an 
approach  negative  answers  may  result,  and  this  situation  may  also  be  handled 
by  linear  programming.  A  condition  v/arranting  the  application  of  linear  pro- 
gramming is  that  there  exist  a  large  number  of  alternative  choices  in  which 
any  one  choice  of  a  subordinate  phase  will  affect  the  outcome  of  the  entire 
plan.  If  there  is  adequate  machine  time  available  in  a  plant  for  a  particu- 
lar product,  there  is  still  an  advantage  in  using  the  method.  However,  if 
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a  number  of  products  can  be  processed  on  any  number  of  machines  in  a 
plant,  a  trial  and  error  procedure  is  far  too  time  consuming,  and  even 
then  the  solution  may  not  be  the  best  one.     The  word  scheduling  is  pur- 
posely omitted  here  for  if  the  problem  involves  a  decision  on  sequence 
of  operations,   present  methods  of  linear  programming  cannot  handle  this 
completely.     Research  is  being  conducted  along  this  line  such  as  noted 
in  reference  (5). 

C names6  has  described  a  linear  programming  problem  as  follows: 
If  a  problem  can  be  written  in  the  following  form, 

Given  a  criterion  of  choice 

n 

(1)  2l  c-jX-<  which  is  to  be  minimized  or  maximized,  and  is  subject 

3=1    rJ 

to  the  tiro  types  of  constraints  or  conditions  that  follow: 

n 

(2)  21    a^^bj  ,  i  =  1,   2,   3....m,  and 
j«l 

(3)  x.  >0  ,  j  =  1,  2,  3...«n;  then  the  problem  is  a 
linear  programming  problem. 

Here,  the  objective  is  the  criterion  of  choice  (1),  and  the  con- 
straints are  contained  in  (2)  and  (3)»  For  maximizing,  the  following 
meanings  may  be  associated  with  the  variables: 

Xj=a  variable  to  be  determined  such  as  the  quantity  of  a  part  or 

a  product. 
c*  =  profit  per  unit, 
a^  r  a  capacity  constraint  such  as  a  unit  production  time, 
b^  =  a  stipulation  such  as  production  capacity  expressed  in  time 
available. 
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Although  solution  to  the  following  problem  is  obvious,   it  is  pre- 
sented to  show  formulation  to  the  above  conditions.     Suppose  there  are 
two  products     P,     and     P     ,  one  machine  available  that  can  produce  ei- 
ther product,  and  the  available  machine  time  is  8  hours. 

x-,  =  units  to  produce  ^2  e  units  to  produce 

2  hrs  =  production  time  per  unit  h  hrs  s  production  time  per  unit 

$l|.00  =  unit  profit  $2,00  =  unit  profit 

Problem: 

(1)  $li#00  x1+  $2,00  x_  ,  maximize  profit 

(2)  2  hrs  x    +  3  hrs  x  <    8  hrs 
x.,  ^.  0  ,  and    x?  ^  0 

Aside  from  product  demands  or  promised  delivery,   the  answer  is  apparent 
by  inspection.     However,  as  the  number  of  products,  machines,   or  produc- 
tion paths  available  increase  and  product  demands  are  considered,  the 
problem  becomes  proportionally  complex.     The  linear  programming  method 
is  carried  out  by  successive  steps,  and  if  a  solution  is  possible,  a 
maximum  profit  case  will  be  reached.     Also,  lower  profit  cases  can  be 
considered  for  factors  such  as  balancing  machine  time  or  other  economic 
reasons. 

Regarding  the  construction  of  the  mathematical  model,  Bellman'   has 
described  operations  of  a  multi-stage  process  to  be  of  two  types:     de- 
terministic and  stochastic.     Deterministic  may  be  considered  as  known 
data,  the  solution  being  completely  determined  for  each  choice  of  al- 
ternative, whereas  stochastic  data  involves  a  probability  distribution 
of  some  form.     References   (27)  and  (28)  discuss  methods  where  probability 
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considerations  are  given  both  to  the  functional  and  to  the  constraints 
involved.  Drawing  an  analogy  in  the  field  of  mechanics  concerning  the 
use  of  statics  or  dynamics,  there  is  as  large  a  field  of  application  in 
Industrial  Engineering  for  use  of  the  deterministic  concept  as  of  the 
stochastic  concept. 

Many  examples  of  linear  programming  emphasize  a  solution  for  the 
number  of  products  to  produce,  ship,  buy,  or  similar  functions.  How- 
ever, one  of  the  most  useful  attributes  of  a  linear  programming  solu- 
tion is  that  of  cost  analysis  and  the  effect  of  changing  constraints  or 
stipulations  such  as  minimum  or  maximum  requirements,  unit  time,  unit 
cost,  unit  rates,  and  similar  factors. 


MATHEMATICS  OF  LINEAR  PROGRAMMING 
This  section  is  a  summary  of  the  basic  mathematical  tools  and 
concepts  involved  in  linear  programming.     It  is  not  all  inclusive  but 
contains  the  essential  features  normally  encountered  in  the  literature 
concerning  this  field.     References  should  be  consulted  for  precise 
definitions  and  details. 

MacDuffe     points  out  that  the  importance  of  matrix  theory  derives 
from  the  rules  of  combination  of  matrices  and  the  fact  that  they  may  be 
represented  as  square  arrays  is  incidental.     The  term  "tableau,"  as  used 
in  linear  programming  for  the  tabular  array,   indicates  the  array  to  be 
a  member  of  matric  algebra.     The  tableau  is  also  referred  to  as  a  frame, 
since  it  describes  the  convex  hull  of  the  polyhedral  cone  by  means  of 
half  lines  spanning  the  cone. 

Notation 
Summations . 

A  series  of  values,  known  or  unknown,  may  be  given  as: 

alxl+  a2x2  +  a3X3  ■*"  %Xh  +  a5X5   "•■*"  an*n 

Or  more  compactly: 

n 

2L  a.x,  j  =  1,  2,  3,   ...n 

j=l     J  J 
For  a  set  of  expressions: 

a11x1+  a-^g   ...  4-  a^.  ...  +■  a^  =    £  a^  i  =  1 

J  — X 

n 

a21xl+  ^2  •••  f  a2jxj   •••  f  Vn  =    .^  Vj    -  i  =  2 


ailxl*  ai2x2   —  +  alfl   •••  +■  ainxn  =   ^  aij*J  *  =  rowv 

°  °  j=l       d  °  number 

amlxl+am2x2   —  +  VCJ   —  ^  amnxn  =    j£  aifj  i  =  m 


The  matrix  of  this  set  may  be  denoted  by    [a     ]     or     A   ,  with  each  ele- 

■*■  J 

ment  a.  .  corresponding  to  the  element  in  the  1th  row  and  J*"11  col- 


If  the  column  values  of  each  x.  are  not  identical,  the  above  may 

J 

be  extended  to  use  of  a  double  summation: 

m   n  Where  the  corresponding  elements,  a    and  x.  .  , 

21  21  a.  "*■  ij       iJ 


j*u 


i=l  j=l  °     °        are  multiplied  and  summed  for  total  value, 
In  matrix  notation  the  set  of  expressions  may  be  shown  as: 
n 


U..3 


js-1  ^  J 


i  =  1,  2,  3- 


Inequalities  and  Equalities, 

Symbols  used: 

y  >  x  y  is  greater  than  x 

y  ^  x  y  is  greater  than  or  equal  to  x 

x  <  y  x  is  less  than  y 

x  ^  y  x  is  less  than  or  equal  to  y 

These  symbols  represent  such  restrictions  as  "produce  no  more  than," 

"stock  no  less  than,"  "pay  no  more  than,"  "produce  at  least,"  etc.  In 

the  above  examples  the  latter  two  were  obtained  by  reversing  the  symbol 

and  order  of  variables.  However,  in  solving  a  problem  containing  a  set 

of  mixed  inequalities,  it  is  necessary  to  make  the  system  consistent. 

This  is  carried  out  by  multiplying  through  by  minus  one,  which  reverses 

the  inequality  and  preserves  the  order  of  variables. 

Example:  x    +  x    +■  x _  >y 

Xl+  X2+  X3^  Z 

Depending  on  the  problem  being  solved  (maximize  or  minimize),  either 
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inequality  may  be  multiplied  through  by  minus  one,   giving: 


xl  "  xo  "*  xt  ^ 


xl  "**  x2  "*"  X3  ^  z 


x,  4-  x0  •+•  x-.  ^.y 


x3^ 


Matrices 
Essentially,  matrices  are  sets  of  numbers,  real  or  complex,  ar- 
ranged in  rectangular  arrays  and  which  obey  certain  rules  of  multipli- 
cation. 

Matrix.  A  rectangular  array  of  elements  arranged  in  m  rows  and  n 
columns.  The  matrix  is  of  order  (m,  n).  It  is  written  as  a  (mxn) 
matrix.     If     m  =  n  ,  it  is  termed  a  square  matrix  of  order     n. 


Example : 


A  = 


all 

a21 
a31 

al2 
a22 
a32 

a13  •-  al» 
a23  -  a2n 
a33  •"  a3n 

■     !>«> 

aml 

am2 

a  ,  ...  a„, 

m_>             ran 

Wher 

i  r  1,    2   . 
J  -  1.   2   . 


Where     a,  .     is  each  element, 


Diagonal  Matrix.     A  square  matrix  having  all  off  diagonal  elements  equal 
to  zero.     The  diagonal  elements  may  equal  any  value  except  zero. 
Example : 


1 

0 

0 

0 

3 

0 

0 

0 

2 

Unit  Matrix.     A  diagonal  matrix  having  diagonal  elements  equal  to  unity 
and  off  diagonal  elements  equal  to  zero. 

This  matrix  is  also  symmetric  because 
the  matrix  is  equal  to  its  transpose. 


Example : 

10     0 

I  * 

0     10 

0     0     1 
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Transpose  of  a  Matrix.     The  transpose  of  a  matrix  is  obtained  by  inter- 
changing rows  and  columns  and  is  denoted  by  a  prime. 


Example : 

12     1 

A  = 

3     2     1 

2     3     3 

,     Transpose    s 


1  3     2 

2  2     3 
113 


o  A1 


Row  Vector .     A  set  of     n     numbers  in  a  row,  a  matrix  of  order     (1,  n) 


M 


tai 


'3  ••••an1 


Column  Vector.  A  set  of  m  numbers  in  a  column,  a  matrix  of  order 
(m,  1). 

[b]  = 


Addition.  To  add,  both  matrices  must  have  the  same  number  of  rows  and 
columns,  corresponding  elements  are  added. 

Subtraction.  The  same  rule  applies  as  for  addition  with  corresponding 
elements  subtracted. 

Multiplication  of  Two  Matrices.  Matrix  multiplication  is  defined  such 
that  linear  transformations  may  be  concisely  and  conveniently  expressed. 
Two  matrices  can  be  multiplied  only  if  the  number  of  rows  in  the  first 
matrix  is  equal  to  the  number  of  columns  in  the  second  matrix.  If  this 
condition  is  satisfied  the  matrices  are  said  to  be  conformable.  The  re- 
sult of  this  multiplication  is  a  matrix  with  the  number  of  rows  equal 
to  the  number  of  rows  in  the  first  matrix,  and  the  number  of  columns 
equal  to  the  number  of  columns  in  the  second  matrix.  Elements  of  the 
product  matrix  are  obtained  by  multiplying  a  row  of  the  first  matrix 
by  the  same  numbered  column  of  the  second  matrix. 
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all 

S12 

a13 

X 

321 

322 

a23 

7 

V 

J32 

a 
33. 

z 

[x  y  2] 

"all 

a12     a13 

a21 

a22 

a 

23 

V    11  12y         13   J 

(a21X+a22y  +  a23Z) 
(a31X+a32y+a33Z) 


L  31      32    *33 


=  [(xan-i-  ya21  +  za31)(xan^+  ya. 


za32)(xa13 


'!2-  *~22  -  -«32/v-»13  t-  ya23  +-  za33 
Multiplication  of  a  row  vector  by  a  column  vector  results  in  a  scalar, 


)] 


Cax    a2    a3] 


a  (a-,^  +  a9b5  +  a^) 


klwl 


2W2 


Multiplication  of  a  Matrix  by_  a  Scalar,     Each  element  of  the  matrix  is 
multiplied  by  the  scalir. 

Division  by  a  Matrix.     Division  by  a  matrix  is  performed  by  multiplica- 
tion of  its  inverse. 
Example:       A  B  =  C 


B  =  A"*1  C   ,     and 


=  B 


-1 


Inverse  of  a  Matrix.  The  inverse  of  a  matrix  A  is  written  as  A 


-1 


The  product  of  a  matrix  and  its  inverse  is  the  unit  matrix     I, 

A  A-1  *  A"1  A  =  I 
Elimination  Method  of  Finding  the  Inverse.       For  a  matrix  with  a  number 
of  zeros  present,  the  elimination  method  may  be  simpler  than  the  cofac- 
tor  method.     In  general,   it  should  not  be  used  if  there  are  two  or  more 
identical  columns  present  in  the  original  matrix. 
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Example 


.6 


A  =     Oij} 


A^k)     denotes  the     kth    column  of  matrix     A 
written  in  terms  of     i'*'     and  the  co 
.00 


l(k) 


A  = 


10     0 

2  10 

3  U     1 


efficients  of  V 

denotes  the  k^"  column  of  the  inverse 

matrix  A*1  . 

A(l)  =  i(D-»-2l(2)  +  3l(3) 

A(2)  =       i(2)+.Ul(3) 

A<3>  =  l(3) 


Solving  for  the  I's  in  terms  of  the  A's, 

j(l)  _  A(l)  _  ^(2)  +  5A(3)     Therefore,  A"1 
l(2)a        A(2)-UA(3) 
l(3)  =  A<3> 


10     0 
■2     10 

5  -k    l 


Vectors  and  Points  in  Space 
The  method  of  linear  programming  consists  of  examining  various 
points  in  space.     A  point  in  space  may  be  considered  by  its  coordinates 
alone  or  by  a  position  vector  terminating  at  that  point,  the  position 
vector  being  determined  by  the  coordinates  and  the  origin  of  the  system. 
In  a  3  coordinate  system  (x,  y,  z)  the  unit  vectors  along  each  axis  have 
a  magnitude  of  unity.     These  vectors  may  be  referred  to  as     i,   j,  and  k, 
and  a  point  in  space  may  be  expressed  in  terms  of  these  vectors. 


T" 


Fig.  1    Unit  Vectors 


Ill 

Disregarding  unit  vectors  for  the  moment,  the  segment  joining  two  points 
in  space  in  terms  of  other  vectors  is  shown  in  Figure  2. 


Fig.  2  Line  Segment 


By  analytic  geometry,  the  line  segment  Pg  to  P^  is  described  and 
generated  in  terms  of  vectors  A  and  B  by  mB  +  (l-m)A  ,  where 

Points  in  Space. 

A  point  in  space  may  be  expressed  in  terms  of  the  unit  vectors 

along  each  axis  and  a  scalar, 
z 


Fig.  3  Points  in  Space 


IS 


Let     PA  r  f(x,y,z) 


i,  y  =  2, 


Then  PA  r  xAi  +  yAj  -I-  zAk  =  i  +  2 j  +•  3k 

Let     PB  =  f(x,y,z)  x  =  2,  y  b  U,  z  =  2 

Then  Pg  =  xBi  -I-  yBj  -I-  Zgk  =  2i  +  U j  +■  2k 

PA+  PB  =  (XA+  V1  +  (^A  +  3^   +  (ZA+  ZB)k 

Pc  =  PA+  PB  =  (1+  2)H-(2+-  U)j+-(3+  2)k  =  3i+  6j+  5k 

The  effect  of  multiplying  each  vector  by  a  scalar  is  a  multiplica- 
tion of  each  unit  vector. 

Thus:  3PA  =  3(xAi  +  yAJ  +  zAk)  =  3xAi  +  3yAJ  +■  3zAk  r  3i  +  6 j  -»-  9k 

2PB  =  2(xBi  +  yBJ  +  Zgk)  =  2xBi  +■  2yBi  +  2Zgk  =  Ui  +  8 j  +  Uk 
3PA+-  2PB  r  (3xA+    2xB)i  -»-  (3yA  +  2yB)j   +  (3zA  -I-  2zB)k  = 

=  (3+  U)i  +  (6  +■  8)j   -f-(9-f-  U)k  =  7i+  lUj  +  13k 
Vectors  or  points     PA  ,     Pg    and     Pc     may  be  expressed  as  column  or  row 


or      CxA,  yA,  zA,l     =     [12  3] 
or       [_xB,  yB,  zB,]     r     [2  U  2] 

or      l*C,  7C  ZC>]      -     [3  6  5] 


Because  of  the  large  number  of  variables  involved  in  a  linear  pro- 
gramming problem  and  the  use  of  (1,   2,  3  ...m)  dimensions,  the  axes  are 
expressed  by  numerical  subscripts.     For  a  three  dimension  system  the 
axes  are  normally  expressed  as     x-,    ,     x«  and  x^   ,  for  an     m    dimension 


vectors: 

PA  = 

XA 
*A 

" 

1 

2 
3 

PB  = 

V 

^B 

'2 
h 
2 

pc  = 

"xc" 

ZC 

- 

'3" 
6 
5 
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system  (x,,  x2,  x~,  ....x  ).  The  concept  of  unit  vectors  is  implicit 
and  they  are  not  normally  written. 

In  the  previous  example  the  unit  vectors  may  be  expressed  as  row 


or  column  vectors  in  the  form 


or       [i  j  k] 


i  = 


3  = 


PA  =  i  +  2j  +  3k 


By  matrix  multiplication  and  addition! 


PA  = 


"l" 

0 

0 

0 

-»-    2 

1 

■+-  3 

0 

0 

0 

1 

*-    J 

L  J 

Similarly,   PB  =  2 


1 

0 

V 

1 

0 

■*- 

2 

4- 

0 

= 

2 

0 

0 

3 

3 

l1 

0 

0 

2 

0 

+  h 

1 

■+-   2 

0 

= 

h 

0 

_0_ 

1 

2 

V 

"o" 

"o" 

"3" 

0 

4-  6 

1 

+   $ 

0 

= 

6 

0 

0 

1 

5 

Also, 


PC  =  PA+PB  = 


12  3 

2      +      U       =      6 
_3j  2]  [* 


Chames-*-  denotes  two  spaces  used  in  linear  programming,  a  "re- 
quirement space"  containing  the  required  condition  vectors  and  the 


17 
"solution  space"  containing  the  solution  vectors.     By  means  of  linear 
transformation,   it  is  possible  to  pass  from  one  space  to  the  other  pro- 
vided the  transformation  of  unit  vectors  is  known  or  equivalently,  the 
transformation  matrix.     If     x^,  x2,   •  ••*n    are  the  original  variables, 

and    x',  x»,    ...x1     are  the  new  ones,   they  may  be  exraressed: 
l      c  n 


Xl  =  allxl  +  a12x2 


xn  =  anlxl  +  an2*2   ••••  annxn 
and  the  transformation  matrix  is: 


alla12  ••••  aln 

_anlan2  ••*•  ann_ 
Denoting  the  transformation  matrix  by    A  ,   it  may  be  written: 


x» 

xl 

=  A 

V 

X1 

X2 

X2 

• 
• 
• 

X' 

n 

• 
• 

and  solved  for  in  terms  of  either  set  of  variables, 


zA-1 


Letting  A"1  be  called  a  new  matrix  B  with  elements  b 


ij  ' 


x 

1 

X2 

= 

b    x'  + 

11  1 

b21Xi  + 

b    x» 

12  2 

b22*2 

•  •  •  • 
. . . . 

+-  b    x' 
In  n 

X 

n 

= 

b  ,x»  ■+- 
nl  1 

bn2*2 

• » . . 

+  V<A 
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Convex  Sets 

Theorems  and  proofs  regarding  convex  sets,  convex  cones,  and  convex 
polyhedral  cones  are  given  by  Charnes,1  Blackwell  and  Girshick,^  Gale, 
and  Gerstanhaber.-^-     Gale  states  that  the  theory  of  polyhedral  cones  is 
equivalent  to  the  theory  of  inequalities  with  the  former  being  a  geo- 
metric interpretation  of  the  latter, 

A  convex  set  of  points  is  such  that  if  two  points  are  in  the  set 
then  the  line  segment  joining  them  is  also  in  the  set.     Examples  of  con- 
vex sets  are  points,   lines,  spheres,  hyperplanes,  and  half  spaces. ?     Hy- 
perplanes  are  described  by  equation  (2),  and  half  spaces  by  equation  (3)« 

Theorem:       A  linear  transformation     L    from  an     n-     dimensional 
space    W     takes  a  convex  polyhedron     K     into  a  convex  polyhedron 
L(K)   ,  the  image  of     K. 

A  complete  treatment  of  this  transformation  is  given  in  reference 

(l)j  the  following  figure  is  taken  from  that  publication. 


U  (n-dim.) 
Solution  Space 


W  (m-dim.) 
Requirement  Space 

P-, 


Fig.  h     Solution  and  Requirement  Space 


L(K) 


The  configuration  of  L(K)  is  determined  by  the  transformation 
matrix  but  will  always  be  an  image  of  K  .  Here  the  convex  set  K  of 
3  dimensional  space  is  transformed  into  L(K)  in  2  dimensional  space, 
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the  points  e-,,  e2,  and  e^  being  transformed  into  points  P,,  P„,  and 

p3. 

A  convex  polyhedral  cone  is  generated  by  the  half  lines  out  of  the 
origin,  or  rays,  and  the  convex  set.  Fig.  3  contains  a  convex  polyhe- 
dral cone  with  vertex  at  the  origin  and  the  triangular  convex  set  formed 
by  points  P^,  Pg,  and  Pq  .  Koopmans   gives  two  equivalent  definitions 
of  convex  polyhedral  cones: 

Such  a  cone  Cconvex  polyhedral]  can  be  defined  either  as  the 
convex  hull  of  a  finite  number  of  ha If -lines  out  of  the  origin 
(sum  definition)  or  as  the  intersection  of  a  finite  number  of 
halfspaces  whose  bounding  hyperplanes  pass  through  the  origin 
(intersection  definition). 

An  extreme  point  of  a  convex  set  is  a  point  which  does  not  lie  in 
a  segment  joining  two  other  points}  they  are  the  corners  of  a  set  such 
as  (e-p  eg,  eJ  or  (P-^,  P2,  P?)  of  Fig.  U.  The  extreme  points  are  the 
points  at  which  the  linear  functional,  the  objective,  takes  on  its  maxi- 
mum (or  minimum)  value.  The  general  form  of  linear  programming,  the 
simplex  method,  provides  a  systematic  method  of  examining  every  extreme 
point  and  of  obtaining  the  solution  required  by  the  objective. 

Basis 
A  vector  space  of  m  dimensions  is  spanned  or  generated  by  a  set 
of  linearly  independent  vectors.  The  number  of  linearly  independent 
points  or  vectors  for  an  m  dimensional  space  is  at  most  m  and  this 
set  of  points  or  vectors  constitutes  a  basis  of  that  space.  A  basis 
must  consist  of  a  set  of  vectors  such  that  every  vector  in  the  space 
can  be  written  as  a  unique  linear  combination  of  the  set.  Thus  in  Fig. 
2,  the  space  generated  by  the  rectangular  coordinate  unit  vectors  has  a 
basis  of  its  unit  vectors.  Also,  vectors  A  and  B  can  be  considered 
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as  a  basis  for  the  grid  area  bounded  by  vectors     A     and     B   ,  as  all 
points  on  that  grid  space  can  be  uniquely  described  in  terns  of  the  two 
vectors  and  scalars. 

The  basis  provides  an  initial  starting  point  in  the  solution  of 
linear  programming  by  the  simplex  method.     The  number  of  linearly  inde- 
pendent vectors,   the  rank  of  the  matrix,  determines  the  number  of  di- 
mensions involved  in  the  requirement  space.     In  practice  the  number  of 
dimensions  of  the  basis  in  the  "simplex  method"  is  determined  by  the 
number  of  restrictions  that  must  be  met.     For  example,   in  equation  (2), 
this  would  be  determined  by    m  •     In  the  case  of  the  transportation 
method  the  number  of  variables  in  the  basis  is     n  +  n-  1  . 

Slack  Vectors 
The  general  form  for  a  linear  program  problem  was  given  with  the 
objective  and  restrictions  expressed  in  terms  of  inequalities.     The 
system  of  inequalities  is  changed  to  a  system  of  equalities  by  adding 
new  non-negative  slack  variables,   sometimes  called  giveaway  factors  or 
disposal  activities.     In  matrix  form,  the  slack  variables  form  slack 
vectors,  and  yield  an  identity  or  unit  matrix  which  provides  an  initial 
basis.     Similarly,   if  the  objective  and  restraints  are  set  up  in  the 
form  of  equalities,  artificial  vectors  may  be  used  to  provide  a  basis. 
In  the  case  of  a  maximizing  objective,  a  large  negative  "profit"  value 
may  be  associated  with  the  artificial  vectors;  the  value  itself  need 
not  be  specified. 


GENERAL  FORM  OF  LINEAR  PROGRAMMING 

Simplex  Method 
The  general  forn  of  solution  for  any  linear  programming  problem  as 
developed  by  Dantzig,^  is  the  simplex  method,     Charnes  and  Lemke  •*  de- 
veloped the  modified  method,  a  procedure  better  adapted  to  machine  com- 
putation.    Special  cases  of  the  general  method  have  been  developed  such 
as  the  transportation  method  and  the  assignment  method. 

The  general  form  as  stated  was: 

n 

(1)  Maximize  21    c-^i 

0=1    JJ 

n 

(2)  Subject  to      ZL  a.  jc  .  ^  b.  i  =  1,  2,  3  ••••m 

(3)  Xj  >0  j  =  1,  2,  3  ....n 

For  m  =  3     and  n  =  3  »  equations  (1)  and  (2)  take  the  following  form: 

(h)  Maximize  cixi  "*"  C2X2  "*"  C3X3  =  z 

(5)  Subject  to  a^ixi  +  ai2x2  "*"  a13x3  ^  bl 

a31xl  +  a32x2  *  a33X3  *  b3 

(6)  xj  ^0  j  r  1,  2,  3 
Adding  slack  variables,  inequalities  of  (5)  become: 

(7)  a^  ♦-  ai2x2  +  a13x3  *  a^xu  +  a-^  +  al6x6  =  t^ 

a21xl+  a22x2  *  a23x3  *"  a2liXU  *  a25x5  +"  a26x6  =  b2 

a31xl  +  a32x2  *  a33x3  *  a3UxU  *  a35x5  +  a36x6  =  b3 

x.  >,0 

Equation  (7)  is  the  form  used  in  the  simplex  method  of  solution.     The 

addition  of  slack  variables  also  changes  the  linear  functional  (h)  to: 
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(8)  Maximize    xici"*"  X2C2  *"  X3C3  +  XUCU  +  x£c5  *~  X6C6  =  z 

Before  proceeding  to  a  numerical  example,  features  of  (7)  and  (8) 
will  be  discussed j  the  terms  used  were  previously  described.  The  prob- 
lem is  to  find  the  values  of  x-,,  Xp,  x-,...*X/-  which  maximize  the  linear 
functional  z  .  There  exists  an  infinite  number  of  solutions,  and 
through  use  of  the  simplex  method  the  extreme  points  of  the  convex  set 
are  examined  systematically  in  such  a  manner  that  the  value  of  z  is 
obtained  for  a  maximum  in  the  final  tableau,  frame,  or  table. 

The  column  coefficients  s;q>  a?l>  a31>  ••••>  ai6*  a26*  an<^  a36 
represent  vectors  in  the  requirement  space  W  and  the  row  values 
x,,  Xp,  x~,  Xi  ,  Xj-,  X/,  represent  a  vector  in  the  solution  space  U. 
The  column  coefficients  ai-  may  be  written  as  P^,  P2  ....  P5  and 
the  constraints  b-^,  b2,  b^  as  P  .  Thus,  (7)  may  be  written  as: 

(9)  PlXl  +  Ppx2  +  P3X3  +  Puxu  +  P5x5  +■  P6x6  =  P0 


With  P, 


a2j 

a... 

3j 


Any  set  of  X-j's  which  satisfy  (9)  is  a  feasible  solution  and  the 
set  maximizing  the  value  of  z  ,  in  (8),  is  the  desired  optimal  solu- 
tion. The  values  of  each  a^i     and  0*  may  be  positive,  negative  or 
equal  to  zero.  For  an  objective  of  maximizing,  meanings  which  may  be 
associated  with  the  a^j,  c*,  and  x^  were  given  previously  for  illus- 
tration, but  they  are  not  restricted  to  that  interpretation.  As  first 
given,  the  objective  was  to  maximize  profit  with  z  in  terms  of  dollars, 
but  the  method  can  be  extended  to  any  unit  of  measurement,  time,  pounds, 
ton  miles,  etc.,  and  involving  purchasing,  production,  stocking,  selling 
or  any  similar  function. 
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The  development  and  proof  of  the  simplex  method  in  reference  (2) 
was  made  under  the  assumption  of  nondegeneracyj  degeneracy  is  the  con- 
dition of  having  fewer  than  m  positive  values  of  the  unknown  variable 
at  some  stage  of  the  iteration  process.  This  restrictive  assumption 
■was  first  resolved  for  the  transportation  method  and  later  resolved  for 
the  general  case  by  Charnes.-*-*  •"4 

The  dual  problem  to  the  general  case  may  be  written: 

m 

(10)  Minimize  ^   wibi 

isl 

m 

(11)  Subject  to        51  •*iai.t  >c*  j  =  1,  2,  3  ....n 

i=l  J         J 

w±  >0  i  =  1,   2,  3  ....m 

Here,  the     a±j»  *>±  a^d  Cj     have  the  same  meaning  as  in  the  original 
problem  and  the    w^    would  be  in  dollars  per  unit  time,  a  production 
cost.     Note  the  change  in  direction  of  inequalities  and  the  change  of 
position  for  the     b^     and     c=   •     Here,  slack  vectors  with  a  high  pen- 
alty cost  would  be  introduced  in  lieu  of  the  original  procedure.     The 
general  method  of  solution  is  the  same. 

Example  of  the  Simplex  Method 
A  problem  that  well  illustrates  the  versatility  of  the  simplex 
method  is  the  production  problem  of  references  (15)  and  (26),     The 
problem  consists  of  producing  two  products     R     and    S     by  a  two  stage 
process,   initial  operations  must  be  performed  in  Machine  Center  I,  and 
final  operations  may  be  performed  in  Machine  Center  IIA,   IIAA  or  IIB. 
Machine  Center  IIA  and  IIB  differ  in  that  different  unit  production 
rates  and  costs  are  involved.      In  addition  an  amount  of  overtime  is 
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available  in  IIA  and  is  handled  separately  as  IIAA.  Unit  production 
time,  unit  profit,  and  machine  hours  available  per  machine  center  are 
given  in  Table  1. 

Table  1.  Production  Data 


Machine 
Center 

Product  R 

Product  S 

Maximum  No. 

Operations 

Rl 

R2 

R3 

Sl| 

s5 

s6 

Hours  Avail. 

Initial 

I 

.01 

.01 

.01 

.03 

.03 

.03 

850 

Final 

IIA 

IIAA 

IIB 

.02 

.02 

.03 

.05 

.05 

.08 

700 
100 
900 

Unit  Profit 

.Uo 

.28 

.32 

.72 

.6h 

.60 

There  are  3  ways  of  producing  R,  and  3  "ways  of  producing  S;   therefore, 
each  way  is  considered  to  be  a  different  product  for  computation  pur- 
poses. Consider  the  initial  problem  to  be:  maximize  the  total  profit 
utilizing  the  available  data.  There  are  no  minimum  quantities  of  either 
product  at  this  time;  a  requirement  of  this  type  would  be  handled  as  an 
additional  constraint.  Let  x^,  X2,  x^,  xl,  xt,   and  x^  be  the  unknown 
quantities  to  produce  of  products  R-^,  R2>  R^>  Sl,  St,   and  S^  respec- 
tively. The  problem  may  be  stated: 

Find        xx,  x2,  Xy   x^,  x^,  x6, 


(12)     Maximize 

Z    c1xj  =  z 

(13)     Subject  to 

£aijXj<b 

x^0 

i  =  1,  2,  3 
j  =  1,  2,  3,  U,  5,  6 
The     a^'s     are  the  unit  production  times  of  Table  1,  the  b^'s     are  the 
available  machine  hours,  and  the     Cj's     are  the  unit  profit  for  each 
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product.     Algebraically  the  problem  is: 

Maximize  •hOx1  -f-  .28x2  4-  .32x^  +  «72x^  +  .6bc£  +■  .60x6  r  z 

Subject  to 
(15)  .01x1  4-  .01x2  4-  ,O0x3  +  .03x1|  +•  .03x^  +  .03x6  ^  8£0 

(IS)  .02xx  "♦"•Q&lt  ^700 

(17)  .02x2  4-.0$X£  ^  100 

(18)  .03x3  +.08x6^900 
The  last  four  inequalities  are  changed  to  equalities  by  adding 

slack  variables,     Xj  to  (1$) ,  xg  to  (16),  xQ  to  (17),  and  x10  to  (^K 
The  slack  variables  may  be  considered  as  fictitious  products  with  cor- 
responding    c-'s   ,  unit  profit,   of  zero  dollars  and  therefore  will  be 
forced  out  of  the  desired  optimal  solution  to  (lh).     Equations  (1>) 
through  (18)  may  be  rewritten  as: 

(19)  PlXl  +■  P^  +-  P3X3  +  P^xu  +-  P5x5  4-  P6x6  +-  ?^7  +■  P8x8  4 

4     ?^?  -t-  P^xo  =  P0 

Xj  ^  0  j  =  1,   2   ....10 

10 
or  ^l    PjXj  =  P0         The     Pj's     are  column  coefficients  of  the     a±^ 

matrix  of  (13)  and  of  the  slack  variables. 

The  data  may  now  be  put  in  the  form  of  the  simplex  tableau,  Table 
2.     P0     has  been  placed  on  the  left  hand  side  of  equation  (19)  for  com- 
putational purposes.     Equations  (15)  to  (18),  or  equivalently  (19),  are 
in  the  form  of  (7),  where  some  of  the  coefficients     a^j     are  equal  to 
zero.     The  unit  matrix,  formed  by     Py,  Pg,  PQ,  and  P^q,     is  the  unit 
basis  and  provides  a  starting  point  in  the  simplex  method.     Here  the 
requirement  space  is  four-dimensional,  and  the  solution  space  is  ten- 
dimensional.     Henceforth,  zeros  in  a  tableau  are  indicated  by  a  blank 
space. 
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Table  3  is  formed  from  the  preceding  tableau,  Table  2,  by  repeat- 
ing it  and  listing  the  basis  vectors  in  the  Basis  column,  with  the  cor- 
responding   Cj     for  each  vector.     The  quantities  of  each  product  to 
produce  are  listed  in  the     P0     column.     Thus,   the  basis  automatically 
provides  an  initial  solution;  make  850,  700,   100,  and  900  units  of  the 
fictitious  products     x<7,  xg,  x0,  and  x^q    and  none  of     x^,  X2,  Xo,  and 
xj^   ,  resulting  in  a  net  profit  of  zero*     In  effect,  one  extreme  point 
of  the  convex  set  has  been  examined ;  the  next  step  is  to  examine  the 
other  extreme  points  until  the  objective  has  been  obtained. 

Two  additional  row  values  are  inserted  in  Table  3-     Zj  and  Zj-Cj. 
Extreme  points  are  examined  by  inserting  various  vectors  into  the  basis 
and  testing  each  tableau,  by  means  of  the     Zj- ca     value,   for  an  optimal 
solution.     The  value  of  each     Cj     is   given  and  its  corresponding     z* 
is  computed  from    Zj  =  1  c^ j  ,  where     c±     is  the  coefficient  of  the 

P's  in  the  Basis  column. 

The  following  conditions  for     Zj-Cj     and  the  algorithm  used  in  the 
simplex  method  are  stated  without  proof j  derivations  are  contained  in 
references   (1)  and  (2). 

Conditions  of  z-i-c.:  : 

I  If  any    z.-c.,  <0  ,  and  for  every     i,  Xi,  ^  0     then     z     is  infinite. 

II  If  any    z j-c *  <  0  ,  and  for  some     i,  Xjj  >  0  ,  the  objective  has 
not  been  obtained  and  the  simplex  procedure  is  continued. 

III  If  all     Zj-Cj   ^0  ,  then  the  value  of     z0     is  a  maximum  and  the  ob- 
jective has  been  obtained. 

Conditions  II  and  III  are  normally  encountered  as  condition  I  may 
be  excluded  when  setting  up  the  problem.     In  continuing  the  simplex 
procedure,  following  condition  II,  each  successive  value  of  the  linear 
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functional     z     increases  until  the  objective  is  obtained.     It  should  be 
noted  that  if  there  exists  a     zj-c-j  =  0  ,  for  some     Pj     not  in  the  basis 
of  the  final  answer,  there  are  additional  combinations  of  real  products 
available  but  whose  objective  value  will  equal  the  present     zQ  .     For 
all     P's     in  the  basis,     z-c  =  0     is  always  true. 

The  next  step  is  to  bring  into  the  basis  some  new  column     P^     vec- 
tor.    Any     Pj     having  a  negative     z-c     is  admissible,  but  the  procedure 
is  to  bring  in  that     P     having  the  largest  negative     z-c     value  such  as 
Pi     •     To  determine  which  row  vector,     Pr   ,  to  remove  from  the  basis,  di- 
vide each  element  in  the     PQ     column  by  each  corresponding  positive     a^j 
appearing  in  the     P^     column  coming  into  the  basis.     The  minimum  positive 
ratio  so  formed  determines  the  row  vector  to  be  removed.     This  ratio  must 
be  positive  since  the  original  conditions  required  the  components  of  the 
PQ     vector  to  be  non-negative,  equations   (19). 

For  illustration  the  following  general  tableau  is  given: 


Vector  being 
removed 


Vector  coming 
in 


Vector 

< 

coming 

Basis 

Po 

PJ 

Pk 

pi 

*i 

aij 

aik 

pr 

xr 

arj 

^k 

pi 

alo 

ah 

aik 

pk 

ark 

ark 

ark  =  1 
^k  " 

Old  Tableau 


New  Tableau 


Let     Pj 


any  column  vector. 


P^     =  column  vector  coming  into  the  basis. 
Pr     =  row  vector  being  removed  from  the  basis. 

a^  =  coefficient  determined  by  the  intersection  of  the     P^     column 
and     Pv.     row. 
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In  setting  up  the  successive  tableaus,  the  vector  coming  into  the 
basis  is  determined  first.  P^  and  Pr  are  determined  as  previously- 
stated  and  a-p^  follows  immediately.  The  Pk  row  may  be  written  di- 
rectly by  dividing  through  the  old  entire  Pr  row  by  a^  . 

All  other  rows  may  be  determined  by  the  following  formula  or  al- 
gorithm. Derivations  are  contained  in  references  (1)  and  (17). 
Let  a! .  r  any  other  new  element. 

a^.  =  corresponding  old  element. 

arj  s  old  element  in  the  same  column  and  in  the  row  removed, 
a^  =  old  element  in  the  same  row  and  in  the  column  entering. 
ark  =  old  element  in  the  row  removed  and  in  the  column  entering. 

aij  =  a«"*rJ  ^    (i  *  k) 

The  new  tableau  is  completed  using  the  above  algorithm,     Zj's     and 

the  corresponding    Zj-Cj     are  computed  and  inserted  in  the  tableau. 

Also  the  new  profit  value  is  given  by    z'  =  z  -  9  (zi.-ck)     where: 

X 

the  minimum  positive  ratio  previously  described 


ik 


in  determining  Pr  . 

The  maximum  negative  value  of  Zj-c.:  determines  the  next  Pk  j 

the  process  is  repeated  until  such  time  as  there  are  no  negative  values 

of  z-pC-j  t   and  the  final  tableau  contains  the  optimum  solution.  Also 

the  new  z  ,-c  .  ,  z'.-c  .  ,  may  be  obtained  from: 
J  J     J  J 

rk 
In  summary  the  z.-c.  are  needed  to  determine  the  P,   to  come  into 
the  basis.  By  choosing  the  most  negative  value  of  Zj-c.  ,  it  often  re- 
sults in  a  quicker  solution,  but  not  always.  Knowing  Pk  ,  it  is  then 
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necessary  to  determine  the     ?r     to  be  removed  from  the  basis.     These 
two  phases  are  the  key  to  almost  all  special  methods  of  linear  pro- 
gramming. 

Continuing  the  example,  the  next  tableau,  Table  U,   is  constructed. 
The  most  negative     z-i*<j     is  column     ?\,   ,  the     P^     indicated  by  the 
vertical  arrow.     The     Pr   ,  row  to  be  removed,   is  determined  as  follows: 

min.  JL_  ,       a^  >  0 
1       aik 

850 


■^  r  28,333  for  p?  ,    1°£  r  H;,000  for  P6 


Therefore,     Pq   ,  the  minimum  ratio  row,  becomes  the     Pr  .     Column  head- 
ings    PQ  ....P10    and  the     c.i     values  in  the  corresponding  columns  re- 
main the  same  and  are  omitted  in  succeeding  tables  for  brevity. 

Each  element  in  the  old     Pg     row  is  divided  through  by  the     a  ^ 
value  of   .05,  yielding  lh,000  in  the     PQ     column,    Jk  in  the     P^     column, 
1  in  the     P^     column,  and  20  in  the     Pg     column. 

The  elements  in  all  remaining  rows  are  computed  by  the  algorithm, 

a!  .     =  a .  . TjJ  a.,     : 

ij  ij       ark     ik 

a'     -  a       -380a 
a70        70  "  agf     7U 

a70  =  850-T7o§7-03  =  U3° 


a71-71-4la7U 


371  =  -01  "  Iff)  -°3  =  *  -002 
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Each  element  in  the  new  table  is  determined  as  above.  It  should  be 
noted  that  use  of  the  algorithm  sets  up  a  square  pattern  with  the 
initial  ark  determined.  The  pattern  is  indicated  as  follows: 


aik 


Essentially  the  algorithm  states:  The  new  element  equals  the  old 
element  minus  the  product  of  the  diagonals  divided  by  the  a^  . 
Element  a^g  in  Table  k   is  obtained  from  Table  3  by  the  above  pat- 
tern: 

a78  =  0-U£3^l)=-.6 

The  diagonal  product  may  be  left  to  right,  right  to  left,  and  up 
or  down  depending  upon  the  relative  position  of  the     ark  •     The  result- 
ing tableau  is  given  in  Table  U.     A  shortcut  to  avoid  computing  each 
element  when  zeros  are  present  will  be  discussed  later. 

The     z^»s     are  computed  for  each  column,  and  the  corresponding 
Cj's     are  subtracted  yielding  the     z-r*^     row.     The     zf<M     ar© 
scanned  for  the  most  negative  value,  and  the  process  is  repeated  un- 
til Table  9  has  been  obtained,  the  optimum  solution.     Thus,  the  solu- 
tion calls  for  the  number  of  products     R]_,  R2,  R3     as  indicated  in  the 
PQ    column  with  the  corresponding  profit  in  the     Zj-Cj     entry.     It  also 
calls  for  150  units  of  product    X?  ,  but  this  was  a  fictitious  product 
and  is  ruled  out.     Product    S     is  excluded  since  it  does  not  appear  in 


V 

.ho 

.28 

.32 

.72 

.61* 

.60 

Basis 

po 

pl 

p2 

P3 

*k 

p5 

p6 

P7 

p8 

P9 

P10 

850 

.01 

.01 

.01 

.03 

.03 

.03 

1 

0 

0 

0 

700 

.02 

0 

0 

.05 

0 

0 

0 

1 

0 

0 

100 

0 

.02 

0 

0 

.05 

0 

0 

0 

1 

0 

900 

0 

0 

.03 

0 

0 

.08 

0 

0 

0 

1 

r 

.1*0 

.28 

.32 

.72 

.61, 

.60 

Basis 

Po 

pl 

P2 

P3 

\ 

p5 

p6 

p7 

p8 

P9 

P10 

0 

p7 

850 

.01 

.01 

.01 

.03 

.03 

.03 

1 

0 

p8 

700 

.02 

.05 

1 

0 

P9 

100 

.02 

.05 

1 

0 

p10 

900 

.03 

.08 

1 

M 

zr=j 

-.Uo 

-.28 

-.32 

-.72 

-6U 

-.60 

p7 

1*30 

-.002 

.01 

.01 

.03 

.03 

1 

-.6 

.72 

\ 

]]»,000 

J» 

1 

20 

P9 

100 

.02 

.05 

1 

p10 

900 

.03 

.08 

l 

zi 

10,080 

.288 

.72 

ll4.ll 

z3T°i 

10,080 

-.112 

-.28 

-.32 

-.61, 

-.60 

Hull 

P7 

370 

-.002 

-.002 

.01 

.03 

1 

-.6 

-.6 

.72 

?u 

H,,000 

.U 

1 

20 

,6U 

p5 

2,000 

.1» 

1 

20 

P10 

900 

.03 

.08 

1 

M 

11,360 

,.288 

.256 

.72 

,6U 

3i.lt 

12.8 

■r"3 

11,360 

-.112 

-.021^ 

-.32 

-.60 

Hi.li 

12.8 

cr 

.Ho 

.28 

.32 

.72 

.61, 

.60 

Basis 

po 

pl 

p2 

P3 

n, 

p5 

p6 

p7 

p8 

P9 

p10 

p7 

32.5 

-.002 

-.002 

-1/800 

1 

-.6 

-.6 

-.375 

.72 

Pl* 

3ii,000 

.u 

i 

20 

.61, 

p5 

2,000 

.li 

1 

20 

.60 

p6 

11,250 

.375 

1 

12.5 

z.1 

18, 110 

.288 

.256 

.72 

.61, 

.60 

31*  .1* 

12.8 

7.5 

ZJ^J 

18,110 

-.112 

_.02U 

-.095 

m.u 

12.8 

7.5 

I 

p7 

102.5 

-.002 

-1/800 

.005 

1 

-.5 

-.6 

-.375 

pk^ 

.Uo 

pl 

35,000 

1 

2.5 

50 

-6It 

p5 

2,000 

Jl 

1 

20 

?r*- 

.60 

p6 

11,250 

.375 

1 

12.5 

z.1 

22,030 

.h 

.256 

.225 

1 

.61* 

.60 

20 

12.8 

7.5 

zr>j 

22,030 

-.02li 

-.095 

.28 

20 

12.8 

7.5 

P7 

jl,o 

-.002 

.005 

1/300 

1 

-.5 

-.6 

-1/3 

.1*0 

pl 

35,000 

1 

2.5 

50 

.61* 

p5 

2,000 

.1* 

1 

20 

.32 

P3 

30,000 

1 

8/3 

100/3 

ZJ 

21*,  880 

.1* 

.256 

.32 

1 

.61, 

.853 

20 

12.8 

10.6 

zr*i 

2li,880 

-.021; 

.28 

.253 

20 

12.8 

10.6 

Table  9  -  Final  Solution 


Basis 

po 

pl 

P2 

P3 

n* 

p5 

p6 

'7 

p8 

P9 

P10 

p7 

150 

1/200 

1/200 

1/300 

1 

-.5 

-.5 

-33 

.1*0 

pl 

35,000 

1 

2.5 

50 

.28 

p2 

5,ooo 

1 

2.5 

50 

.32 

P3 

30,000 

1 

8/3 

100/3 

zi 

25,000 

.1*0 

.28 

.32 

1 

.70 

.853 

20 

Hi 

10.6 

aJ^3 

25,000 

.28 

.06 

.253 

20 

H* 

10.6 

3?. 


the  basis.     The  program  is  summarized  in  Table  10, 

Table  10 
Optimal  Program  -  Number  of  Units 


Product  R 


Product  S 


R1  (Centers  I,   IIA).    .  35,000  units 

R2  (I,   IIAA) £,000  units 

R3  (I,   IIB) 30,000  units 


s5 
s6 


Total  profit 


$2^,000 


Table' 11  summarizes  the  program  in  number  of  hours  using  the  data 
of  Tables  1  and  10. 

Table  11 
Optimal  Program  -  Number  of  Hours 


Operation 

Machine 
Center 

Product 

R 

Product  S 

Hours 
Used 

Hours 
Avail. 

Surplus 

Rl 

R2 

R3 

sh    s5  s6 

Hours 

Initial 

I 

350 

5o 

300 

0   0   0 

700 

850 

150 

Final 

IIA 

700 

0 

700 

700 

0 

x 

IIAA 

100 

0 

100 

100 

0 

IIB 

900 

0 

900 

900 

0 

It  should  be  noted  that  when  a  zero  appears  in  the     Pr     row  of  any 
tableau,  the  entire  column  in  which  it  appears  will  be  identical  in  the 
following  tableau.     Similarly,   if  a  zero  appears  in  the     P^     column, 
the  old  row  repeats  itself  in  entirety.     Referring  to  Table  6,  zeros 
appear  in  row    P^     for  columns     P2,  P3,  Pt,   P5,  Py,  Po,  ?iq,     and  th°se 
columns  are  identically  repeated  in  Table  7«     In  Table  6,  zeros  appear 


3U 
in  the  Pk  column,  Pj_  ,  for  rows  ?$   and  P5  •  Thus,  rows  Pk  and  P^ 
in  Table  6  are  identically  repeated  in  Table  7«  Because  of  the  many- 
zeros  encountered  in  a  typical  linear  programming  solution,  the  algo- 
rithm need  not  be  used  for  every  element  when  the  above  conditions  are 
present.  The  general  form,  the  simplex  method,  has  been  criticized  as 
a  time  consuming  process;  however,  by  properly  utilizing  various  short 
cuts,  the  time  involved  is  significantly  reduced.  A  review  of  Tables 
3  to  9  will  quickly  verify  the  above  statement. 

An  inverse  matrix  is  contained  in  each  tableau  if  the  initial  tab- 
leau contains  an  identity  matrix. 

Theorem:   If,  in  an  initial  tableau  we  have  an  identity 
matrix,  then  in  this  position  in  subsequent  tableaus  will 
be  the  inverse  of  the  current  basis. 

Referring  to  Tables  3  and  h,   the  matrix  formed  by  Py,  Pg,  P0, 

and  P-|_q  column  of  Table  h,   is  the  inverse  matrix  of  the  new  basis, 

P7,  P^,  P^,  P1Q  of  Table  U,  as  given  by  the  P?,  P^,  P?,  and  P1Q 

columns  of  Table  3« 


Table  3 


>     p9>     p10 


P7>  PU 

1  .03  0       0 

0  .05  0       0 

0  0  10 

0  0  0         1 


Table  h 

p7>     p8>     P9»     P10 

1      -.6       0        0 

0       20         0         0 

0         0         10 
0         0         0         1 


10  0  0 
0  10  0 
0  0  10 
0     0     0     1 


Similarly,  by  the  same  theorem,  in  each  succeeding  tableau  an  in- 
verse matrix  will  be  in  the  same  position  as  the  previous  identity  ma- 
trix.    In  Table  5  a  unit  matrix  is  given  by     Pr>,  Pr,   ?t,  and  P^g  • 
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Therefore,   in  the  same  columns  of  Table  6,  the  inverse  of  the  new  basis, 
F-7,   IY,  ?t,  P5     of  Table  6,  will  be  given. 


Table  5 

p7>     *V     P5'     P* 

1 


3/100 

0 

0 
8/100 


Table  6 

V     V     P5>     P 

1 


10  _ 
-3/8" 

0 

0 

25/2 


— 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

Carrying  the  theorem  further  leads  to  the  modified  simplex  method 
developed  by  Charnes  and  Lerake. *     This  also  leads  to  the  basis  for  ev- 
ery "adjacent  exteme  point"  method  to  date  and  most  of  the  special  meth- 
ods of  solution.     In  the  modified  simplex  method,  the  algorithm  need 
only  be  applied  to  that  part  of  the  tableau  under  the  unit  vectors,  and 
by  one  equation  the     x^'s     and  the     Zj-c»     may  be  developed.     The     Pi, 
is  then  determined,     a^     and  JL.    developed,     Pr     determined  and  the 

simplex  algorithm  applied  as  before. 

Referring  to  the  last  theorem,   it  is  possible  to  develop  all  ele- 
ments in  subsequent  tableaus  by  its  use.     Using  Table  U  for  example, 


R^Y  +  P^  +  P^9  -t-  P10*10  =  P0 


or  in  matrix  notation, 

[p7  pu  p9  PjoJ 


■10 
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Let  B  =    [P?  fy  P9  Piq]     >  and    x  = 


Then,     Bx  =  P0  = 


x7 
XU 
x9 

J^IO 

850 n 

700 

100 

900 

Multiplying  through  by  the  inverse  of  B,  B""-1-  ,  the  latter  equation 
becomes: 

B""XBx  =  B"^ 


or 


B^P, 


Therefore,  it  is  possible  to  obtain  the  new  x-j^  values  by  means  of  the 
inverse  of  the  current  basis  matrix  in  the  new  tableau  from  the  P0  col 
umn  vector  in  the  old  tableau.  For  example  the  x^  in  the  PQ  column 
of  Table  h  can  be  obtained  as  follows: 

x?     =  1(850)  -  .6(700)  +  0(100)  +  0(900)  =         U30 
x^    =  0(850)  1-  20(700)  +  0(100)  +  0(900)  =  lli,000 

x9 

x10 
Similarly  new  column  elements  can  be  developed  using  the  old     Pj 

column  and  the  inverse  of  the  current  basis  matrix.     Elements  of     P^ 

in  Table  h  illustrate  this: 

ai;L  =  l(.Ol)  -  .6(.02)  +  0(0)  -t  0(0)  s  -  .002 

a2i  =  0(.0l)  +   20(.02)  +  0(0)  +■  0(0)  «       .U 

a31  =  0(.0l)  +     0(.02)  +■  1(0)  +  0(0)  s         0 

aiil  =  0(.0l)  +    0(.02)  +  0(0)  +  1(0)  =         0 


0(850)  +    0(700)  +  1(100)  +  0(900)  z        100 
0(850)  +    0(700)  -1-0(100)  +  1(900)  =        900 
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It  is  to  be  noted  that  the  inversion  of  a  large  matrix  becomes   in- 
creasingly difficult  as  the  size  of  the  matrix  increases. 

The  usefulness  of  linear  programming  in  analysing  cost  data  for 
changes  in  initial  restrictions  such  as  minimum  product  requirements , 
available  machine  time,  unit  production  rates,  and  unit  profits  is  well 
illustrated  by  this  problem. 

Initially,  there  were  no  minimum  requirements  of  product     S     to 
meet  obligations  such  as  market  demands.     Assume  a  minimum  of  8,000 
units  of  product     S     are  required  in  addition  to  the  initial  conditions. 
The  problem  need  not  be  reworked  as  the  necessary  data  is  contained  in 
Table  9,  the  former  solution. 

The     zj-°i     row  represents  "opportunity  costs."     Thus  the  process 
given  by    P^     of  Table  9  has  the  least  value,   .06,  a  value  that  will 
affect  the  total  profit  the  least  for  any  change  to  the  production  pro- 
gram.    The  next  succeeding  values  for  product    S     are  .253  and  .28,  for 
P5  and  Pi      respectively.     Thus,     ?t    represents  the  best  choice  and 
will  necessitate  reducing  the  amount  of    ^     produced,  represented  by 
?2     i*1  the  basis.     However,  the  initial  condition  on    x^     was  that 
Xi  >  0  . 

Therefore,     2.5S^  <.  £000,  and  for  maximum  benefit  the  inequality 
sign  is  dropped  yielding:     2,5S^  =  5000,  or  S^  r  2000  . 

Since  8,000  units  of     S     are  required,  the  remaining  6,000  units 
must  be     S^   ,  S^     or  some  combination  thereof.     The  next  lowest  value 
of  the     2  ■<-£•*     row  is  .253  of     P/-    ,  requiring  a  reduction  of    Xo     pro- 
duced, represented  by     P^     in  the  basis. 

Therefore,     &£  ^  30,000,  and  this  condition  is  satisfied  for 
S6  =  6000.     In  addition,     X3  is  reduced  by:     (8/3K6OOO)  =  16,000   . 
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Table  12  summarizes  the  above  results.     In  the  first  step  the 
profit  was  reduced  by  (2000)(ft.06)  =  $120,  and  in  the  second  step  by 
(6O00)($.253)  =  $3520. 

Table  12 
Optimal  Program  -  Number  of  Units 


Product  R 

Product  S 

R±  (I,   IIA)    .    .    .   35,000 

su  • 

0 

R2  (I,   IIAA).    .   .             0 

%' 

2000 

R3   (I,   IIB)    .    .   .  12^,000 

s6   • 

6000 

Total  Profit  =  25,000  -  (120  +  1520)  =  $23,360 

The     Zj-Zj     row  of  Table  9  may  also  be  utilized  to  analyze  the  ef- 
fect of  changes  to  the  hours  available  in  each  machine  center.     For 
zj"°j   *   J  ~  It   8,  9,   10,   initially  introduced  as  a  result  of  the  h  slack 
variables    xy,  xg,  x^,  and  x^q  ,  the  "opportunity  costs"  represent  val- 
ues of  additional  profit  for  each  additional  machine  center  unit  time 
of  centers  I,  IIA,   IIAA,  and  IIB  respectively. 

From  Table  9  it  is  more  profitable  to  make  additional  time  avail- 
able in  centers   IIA,   IIAA,  and  IIB  as  indicated  by  the     z-i-Cj     of  20, 
Hi,  and  10.6.     150  hours  of  surplus  time  is  already  available  in  center 
I  and  verified  by    Z--C-  =  0  ,     the  least  profitable  center  to  allocate 
additional  time.     Assuming  the  150  surplus  hours  available  are  to  be 
utilized  but  not  exceeded,  the  upper  limits  for  IIA,   IIAA,  and  IIB  are 
obtained  as  follows: 

Center  IIA         150  +  Hg  (-.5)    £  0,  therefore  Hg  ^  300  hrs. 

Center  IIAA       150  +•  Hq  (-.5)    >  0,  therefore  %«  300  hrs. 

Center  IIB         150  +■  H10(-.33)  £  0,  therefore  H10^  U50  hrs. 
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Allocation  of  additional  hours  to  Center  I.IA,    IIAA,  and  103,  within 
the  above  limits,  will  increase  the  original  total  profit  as  follows: 

New  profit  =  25,000  +■  20HX  +  lUHg  •*-  1Q.6H- 
for     0  <=  Hg  ^  300,     0  «  Hq  ^  300,     0  ^  H10  <  k$0 

Change  to  unit  production  rates  of  Table  1,  as  a  result  of  new  tools 
or  machines,  may  require  a  new  solution  of  the  entire  problem.     However, 
the  use  of  the  final  solution  of  Table  9  in  analyzing  such  a  change  can 
be  shown  in  this  problem.     Assume  that  the  machining  time  per  unit  in 
Center  III  and  IIAA  for  product     R     can  be  reduced  from  .02  to  ,0175 • 
The  most  profitable  methods  are     Rp  Ro,  and  R2     respectively.     Allo- 
cating the  total  time  available  in  Center  IIA,   700  hours,  requires  h00 
hours  in  Center  I.     The  next  most  profitable  method  is  for  product     R-j   , 
and  900  hours  of  Center  IIB  are  allocated.     This  requires  300  hours  in 
Center  I  for  product     R-3  .     The  100  hours  of  Center  IIAA  is  allocated 
for  product     R2     which  in  turn  requires  57  «1  hours  in  Center  I.     Table 
13  summarizes  the  new  program.     Profits  are  computed  from  the  number 
of  parts  used  in  determining  Center  I  hours  and  the  initial  data.     The 
total  profit  does  not  take  into  account  the  improvement  cost  but  does 
provide  a  figure  of  S2, 198*80  as  a  criterion  for  judging  the  worth  of 
such  an  improvement. 

Table  13 
Optimal  Program  for  Machine  Time  of  IIAR  =  .0175 


Operation 

Machine 
Center 

Product  R 

Product  S 

Hours 
Used 

Hours 
Avail. 

Surplus 

R-j        Rp       R^ 

\     S5     S6 

Hours 

Initial 

I 

Uoo  57.1   300 

0       0       0 

757.1 

850 

92.9 

Final 

IIA 

IIAA 

IIB 

700 

100 

900 

0 

0 

0 

700 
100 
900 

700 
100 
900 

0 

0 

0 

Total  Profit:     16,000  +■  1598.80  -t-  9,600  =  $27,198.80 


Changes  in  unit  profits  may  or  may  not  require  a  new  solution  of 
the  problem.  In  either  event  Table  9  may  be  utilized  as  an  initial 
starting  point  by  using  the  elements  of  that  table  and  the  new  unit 
profits.  The  new  z  •  and  z ^-<-  are  computedj  if  all  z^-Cj  are 
positive  the  program  is  still  optimum.  If  any  z.-c.-  is  negative, 
the  simplex  procedure  is  carried  out  until  all  such  values  are  positive. 

Assume  that  after  obtaining  the  program  of  Table  9,   it  is  necessary 
to  reduce  the  sales  price  of  product     R     by  $.0U  to  meet  new  competitive 
demands.     Table  lli  is  obtained  from  Table  9  by  repeating  tlie  elements 
and  inserting  new     c='s   ,  $.0U  lower,  for    R-^,  R2,  and  Ro   •     Under  the 
old  program  of  Tables  9  and  lli,  and  the  new  unit  profit  figures,  the 
profit  is  $22,200,  a  reduction   of  $2,800.    With  the  new  program  of 
Table  15,  the  profit  is  $22,280,  a  reduction  of  $2,720  from  that  of 
the  original  solution. 

Table  15  is  obtained  from  Table  lh  in  three  steps.     The     Pk  , 
P^   ,   is  determined  by  the     z-c     value  of  -«0U,  and     Pr     from  the  mini- 

mum  ratio  of    JL.  •     In  Table  15,     P^    row  is  calculated  by  dividing 
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the  old     ?2     row  bY  the     a^     of  2,5.     New  rows     P-j_  and  P_     are  iden- 
tical as  in  Table  lit  in  view  of  the  zeros  present  in  the     P    and  P.    , 
rows  and  columns  rule.     Thus,   the  "square  pattern"  of  the  algorithm 
need  only  be  applied  to  new  row     P«  .     z  .  and  z,-c.     are  then  computed, 
all  values  of  the  latter  are  positive  and  therefore  an  optimal  program 
exists. 
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Table  lit 


lie-* 

•  36 

•2b 

.28 

.72 

.6k 

.60 

Basis 

Po 

pl 

P2 

P3 

Pb 

P5 

P6 

P7 

P8 

P9 

P10 

P7 

150 

1/200 

1/200 

1/300 

1 

-.5 

-.33 

.36 

Pl 

35,000 

1 

2.5 

50 

.2b 

P2 

5,000 

1 

© 

50 

.28 

P3 

30,000 

1 

8/3 

100/3 

ZJ 

.36 

.2b 

.28 

.90 

.60 

.7b6 

18 

12 

9.3 

zr°j 

.28 

-.oh 

.1U6 

18 

12 

9.3 

Table  15 

P7 

lltO 

-1/500 

1/200 

1/300 

1 

-.5 

-.6 

-.33 

.36 

Pl 

35,000 

1 

2.5 

50 

.6b 

P5 

2,000 

.h 

1 

20 

.28 

P3 

30,000 

1 

8/3 

100/3 

zj 

.36 

.256 

.28 

.9 

.6b 

.7b6 

18 

12.8 

9  1/3 

z^j 

.016 

.18 

.U16 

18 

12.8 

9  1/3 

The  usefulness  of  linear  programming  in  cost  analysis  is  an  impor- 
tant characteristic  that  should  not  be  overlooked.     The  manufacturing 
process  involved,  depending  upon  its  complexity,  may  cause  difficulty 
in  normal  cost  analysis.      In  an  oil  refinery  problem,   it  -was  necessary 
to  produce  a  minimum  quantity  of  a  product  for  goodwill  purposes,  and 
at  little  or  no  profit.     Because  of  the  process  involved,   it  was  found 
that  to  produce  only  the  minimum  quantity,  the  resulting  costs  were 
higher  than  they  were  to  produce  a  quantity  above  the  minimum. 
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Reference  (17)   illustrates  a  problem  in  planning  an  experimental 
program,   solved  by  linear  programming  and  game  theory.     Alternatives 
such  as  installing  the  plant  without  preliminary  experimental  tests, 
with  bench  scale,  and  pilot  plant  tests  were  considered. 


TRANSPORTATION  METHOD6'  13»  17>  1? 
The  transportation  method  is  a  special  case  of  the  general  simplex 
method.  It  was  originally  considered  by  Hitchcock  and  by  Koopmans.  As- 
sume there  are  a-^,  a2  . ...a^  homogeneous  products  at  origins  0-^,  02«. 
..0^  ....0m  .  These  quantities  are  to  be  completely  distributed  to  des- 
tinations D-p  D2  ••••Dj  . ...Dn  in  quantities  b1,  b2,  ••••bi  ••••bn  ■ 

m       n 
Total  distribution  requires  that   21  a-?  =  "2L   b^  •  Let  the  unit  cost 

i=l     j=l  J 

of  transportation  from  0^  to  Dj  be  c^a    and  the  amount  to  be  trans- 
ported be  x^j  .  The  objective  is  to  minimize  the  total  cost  of  trans- 
portation, A  special  case  of  the  above  method  is  known  as  the  "assignment 
method"  where  all  of  the  a^'s  and  b.'s  are  equal  to  unity.  The  trans- 
portation problem  takes  the  following  form: 

n   m 
(1)  Minimize    21  21  CifCij 
j=l  i=l 


i  =  1,  2,  3  ....m 


(2) 

Subject 

to 

n 
3=1 

xij  = 

ai 

(3) 

m 

i=l 
xi1 

xij  = 

b3 

j  =  1,  2,  3 


If  a  problem  can  be  put  in  the  above  form,  it  is  possible  to  use 
the  transportation  method.  Written  out  in  general  form  for  three  ori- 
gins and  four  destinations,  equations  (2)  and  (3)  take  the  following 
form: 
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xll  +  x12  +  x13  *  x3li  =  ax 

x21  -v-  x22  +  x23  4-  x^  =  a2 

x31  +  x32  +"  x33  "*"  X3U  *  a3 
xll  x21  x31  =  bl 

x13  x23  x33      =  b3 

xlii  x2l*  X3U  =  \ 

This  could  also  represent  three  products  and  four  machine  operations. 

The  transportation  problem  always  has  the  above  symmetrical  form,  equal 
number  of  rows  and  diagonals.  In  addition,  there  are  exactly  two  co- 
efficients of  unity  in  each  column  and  all  other  entries  are  zero. 

In  a  production  problem  for  example,  the  coefficients  would  rarely 
be  unity  for  production  rates.  However,  it  is  possible  to  change  the 
coefficients  to  unity  and  obtain  an  approximate  solution  by  two  basic 
operations:" 

1.  Multiplication  of  an  inequation  by  a  positive 
constant. 

2.  A  change  in  scale  for  certain  of  the  variables. 

This  procedure  results  in  an  approximate  solution  due  to  the  round- 
ing errors  introduced.  "Normalizing"  is  the  term  used  to  describe  the 
above  operation  which  results  in  w6rking  with  equivalent  production  units 
in  place  of  the  original  production  units.  A  complete  description  for 
handling  the  "normalizing"  procedure  is  contained  in  reference  (19).  It 
also  contains  numerous  examples  as  well  as  a  general  discussion  of  linear 
programming  and  the  transportation  method  in  particular.  As  in  the  sim- 
plex problem,  artificial  vectors  may  be  used  to  handle  problems  where 
the  total  amount  available  does  not  equal  the  total  amount  required. 
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Fictitious  origins  and  destinations  with  fictitious  costs  permit  an  op- 
timal program;  all  fictitious  routes  are  forced  out  in  final  solution. 
By  the  above  process  a  mathematical  analogy  may  be  made,  where  originally 
there  was  no  apparent  physical  analogy.  The  construction  of  the  mathe- 
matical model  requires  complete  familiarity  with  the  different  methods 
of  linear  programming. 

The  following  problem  illustrates  the  "stepping  stone"  method, 

originally  developed  by  Dantzig  (18)  and  described  by  Charnes  and  Coo- 

6  20 
per.  »     it  can  be  shown  that  the  m  +  n  equations,  (2)  and  (3), 

-]  o 

constitute     m  4-  n  -  1     independent  equations  in     m  X  n     unknowns. 
These  equations  therefore  form  a  basis,  and  the  minimizing  solution 
will  involve  at  most     m  +•  n  -  1    routes  with  positive  shipments. 

The  stepping  stone  method  consists  in  setting  up  an  initial  solu- 
tion by  means  of  the  basis.     The     a^     elements  are     -h  1,  -  1,  or  0  . 
Thus  in  tabular  form  the     xi     is  used  directly  in  forming  the  minimum 

positive  ratio    __L.  .     The     c±±  ,  cost  of  each     ij     route,   is  inserted 
aik 


in  the  small  corner  boxes  of  Table  15  for  ease  in  computation*     Quan- 
tities available  are  listed  in  the  right  hand  column,  and  quantities  re- 
quired in  the  bottom  row.     Initially  an  arbitrary  number  is  assigned  to 
some  row,  normally  the  row  containing  the  most  stepping  stone  positions, 
the  vectors  in  the  basis.     The     K^'s     for  stepping  stone  columns  are 
then  determined  from  equation  (ii): 

(U)  Rjl +■  Kj  =  c-ji  where     Rj_  =  Row  number 

K-  =  Column  number 

Row  and  column  values  can  then  be  computed  for  the  remaining  stepping 

stone  positions  by  equation  (U). 
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Row  and  column  values  for  the  non-stepping  stone  positions,   the 
vectors  not  in  the  basis,  can  be  determined  from    R^'s     and     K^'s 
available.     For  the  non-stepping  stone  positions  the     z^-j     is  obtained 
from  equation  ($)t 

(5)  «ij«Ri+Kj 

From  equation  (5),  the     z±^±i     is  computed  for  each  cell  not  in  the 

basis.     The  non-stepping  stone  cell  with  the  largest  positive  value  of 

z^  ^-c .  .     is  then  brought  into  the  solution,  and  the  stepping  stone  cell 

with  the  minimum  positive    xi       ratio  is  removed.     To  bring  a  new  vec- 

aik 

tor  into  the  basis,   it  must  be  brought  in  by  a  path  via  the  stepping 

stone  positions  by  rook  moves,  as  in  chess.  The  amount  brought  in  is 

added  or  subtracted  alternately  to  the  basis  positions  used.  The  amount 

brought  in  is  subtracted  from  the  positive  cell  and  added  to  the  negative 

cell.  In  the  result  of  a  tie  for   i  ,  remove  the  most  northwesterly 
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value,  and  if  the  tie  occurs  in  the  same  row  remove  the  most  westerly 

value. 

New    R.  and  K.     values  are  determined,  the  non-stepping  stone 

Zji-Cji     determined,  and  the  process  repeated  until  all  of  the     z^.j-Cji 

are  negative  or  equal  to  zero.     For  non-stepping  stone  cells  with 

zii"°ii     equal  to  zero,   it  indicates  an  alternate  route  with  the  total 

cost  equal  to  that  of  the  present  solution. 

Example  of  Transportation  Method 
For  a  problem  with  two  origins,  three  destinations,  and  correspond- 
ing transportation  costs  as  given  below,   the  use  of  fictitious  values 


hi 

is  shovm  for  the  case  of  £a^  ^  2Tb^   .     The  objective  is  to  obtain  the 

lowest  total  transportation  cost  and  fulfill  all  requirements  of  the 

destinations. 

60  units  requfd 


$9Amit 


$3/unit 


100  units 
available 


70  units  req'd. 
*7/un2.t  $U/unit 

*•    Dp     -* 

$10/unit 


$5/unit 


U0  units  req'd, 


80  units 
available 


Algebraically  the  problem  takes  the  following  form: 
Minimize         9xvi  +"  T*!?  "*"  -^13  +  3x^t  +  liXpo  +"  ^x?3 
Subject  to     X-JJL+  x-^  "*"  X33  £100 

x21  +  x22  +  x23    ^  ^° 

+•  x2i  =  60 


•11 


-12 


'22 


■13 


+•  x 


23 


70 


=  Uo 


xij>° 


Since  there  are  10  more  units  available  than  are  required,  a  fictitious 
destination  is  set  up,  and  a  zero  transportation  cost  assigned  in  order 
to  insure  such  a  route  is  utilized  for  the  excess  supply. 

With  the  fictitious  destination,  the  above  equations  take  the 
following  form: 
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9x-n  +•  7x19  +■  lCbc-,0  +■  0xiv  +■  3x91  +  l4X09  +-  £x9o  +"  Ox 


"11  T    ,A12 

xll 


■13     T    ^^ 


21  T  ^22 


■23 


~2F 


+•     x-^  +    x-jj 


100 


xll 


x21  "*"   x22  "*"   x23  "*"  X2F  =     $0 

+•  X21  r    60 


k12 


+  x 


22 


x13 


+  x 


23 


■IF 


4-x 


2F 


=  70 
=  UO 
=     10 


Xij^O 


Due  to  the  characteristics  of  this  type  problem,  the  problem  may- 
be reduced  to  the  tabular  form  of  Table  26  for  computation  purposes. 

Table  16 


°1 

x12 

Jfil 

x13 

X1F 

100 

°2 

x21 

x22 

5 
x23 

0 
X2F 

80 

60 

70 

Uo 

10 

180 

Table  17  is  constructed  from  the  data  of  Table  16.  Since  the  number 
of  vector  elements  in  the  basis  is  five  (m-h  n  -  l),  an  estimated  solu- 
tion may  be  written  down  immediately.  For  a  systematic  approach,  par- 
ticularly with  a  large  matrix,  it  is  convenient  to  start  at  the  northwest 
corner,  work  across  and  down,  balancing  requirements  at  each  step.  This 
procedure  is  used  in  Table  17,  total  cost  is  $llU0. 
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Table  17 
h  5 


■4 

60  + 

Uo- 

10J 

±i 

100 
80 

J>J 

30  + 

Ac 

oj 

10 

60 

70 

llO 

10 

180 

Rules  in  establishing  row  and  column  values,  R^  and  K .  ,  for  the 

J 

stepping  stone  positions  are  summarized  below:6 

1.  Arbitrarily  assign  0  to  a  row. 

2.  Move  eastward  along  the  row  until  you  encounter  the  first 
vector  in  the  basis. 

3.  Assign  to  the  column  the  value  necessary  to  cause  the  row 
and  column  sum  to  equal  the  cell  coefficient.     [Equation  (k)~) 

U.     Duplicate  the  process  for  all  cells  in  the  row  which  are 
included  in  the  basis. 

To  establish  row  and  column  values  for  cells  not  in  the  basis:6 

$.     Froceed  to  determine  any  other  row  number  by  selecting  a 
vector  within  the  basis  for  which  a  column  number  is  avail- 
able. 

6.     Carry  through  this  procedure  until  all  row  and  column  num- 
bers have  been  determined. 

Applying  the  procedure  to  Table  17,  assign  zero  to     R2   ,  this  deter- 
mines the     K2,  Ko,  K^     of  h,  5,  0  respectively  by  equation  (U).     R]_    may 
then  be  found  from     K2     and  equation  (h) ,  equal  to  3«     K-i     is  then  de- 
termined by  the  sane  equation.     The     Z-h's     for  ceLL3  13,   lU,  and  21,  are 
found  by  equation  (5);  the  corresponding    Cij's     are  subtracted,   giving 
values  of  -  2,   -t-3,  and  +•  3  respectively.     Cell  21  with  the  largest 
positive     Zj^^-c^^   ,   +- 3>   is  the  vector  to  be  brought  into  the  basis,  and 

cell  22  with  the  lowest  positive  x>       ratio,   30,   is  to  be  removed  from 

aik 
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the  basis.     Here,  two  cells  resulted  in  a  tie  by  the  4-3  value.     In  a 
matrix,  rows  may  be  interchanged  and  columns  may  be  interchanged.     Inter- 
changing row  1  and  2  would  result  in  cell  21  becoming  the  most  north- 
westerly cell,  applying  the  rule  for  case  of  a  tie,  requires  cell  21  to 
be  brought  in.     Cell  21  is  brought  in  via  the  stepping  stones  by  rook 
moves.     Note,     P2^    may  be  written: 

(30)     (U0)     (60) 
P21  =  P22  -  P^   +  Pn 

The     ajj^'s     are   4-1,  -  1,   or  0,  and  the  plus  and  minus  signs  within 
each  cell  refer  to  each    a^     of  the  path.     To  determine     0  ,  only  the 
P's     with  the    +1  coefficients  are  considered,     P22  and  P-q  .     Liinimum 


e  ,  2L  , 


is  4p  ;  the  amount  to  be  removed  is  30  for  a  reduction  in  cost 


of  $3  X  30  or  C?°0.     Total  cost  is  $llii0  -  $90  =  $1050.     The  new  program 
is   given  in  Table  18. 


Table  18 
1  5 


30  + 

70 

ipj 

J>J 

100 

30- 

id 

Ao 

10  4- 

80 
180 

60 

70 

ko 

10 

New  row  and  column  values  are  determined  and  shown  in  Table  18.     The 
Zji- Cj*     for  non-stepping  stone  positions  show  that  Cell  lLi,  a  differ- 
ence of   +6,  must  be  brought  in  through  stepping  stones     P-q,   P2i>  and 
Pr>i      with  minimum     9    at     Pgi    ;  therefore     Pg.      is  removed.     Reduction 
in  cost  is  $6  X  10  =  $60,  total  cost  of  $1050  -  $60  or  $990.     The 


51 


program  is  given  in  Table  19. 


h      9 


Table  19 
7  11  0 


20 -h 

70 

10j 

oj 

10 

100 

Uo- 

id 

Jl, 

Uo  + 

oj 

80 

60 

70 

Uo 

10 

180 

R> 

l 

0 
-  6 


New  row  and  column  values  are  determined,  given  in  Table  19,  the 
z,.<-  differences  computed,  P-jo  is  brought  in,  P^  removed,  and 
the  cost  reduction  of  $20  is  obtained,  $1  X  20.  The  program  is  given 
in  Table  20  at  a  total  cost  of  $970. 

Table  20 
K.     8      7      10      0 


2} 

u 

70 

20 

10 

100 

80 

-4 

60 

kl 

20 

0] 

60 

70 

Uo 

10 

180 

New  row  and  column  values  are  assigned  and  given  in  Table  20.     The 
z^-pCji     for  all  non-stepping  stone  positions  are  negative  and  an  op- 
timal program  has  been  obtained.     If  any  of  the     z^-Cji     far  the  above 
positions  were  equal  to  zero,  such  a  cell  would  represent  an  alternate 
route  but  at  the  same  total  cost. 


Example  of  Transportation  Method,  Fictitious  Route0 
To  illustrate  other  considerations  concerning  the  transportation 
method,  the  following  example  is  given.     Procedural  explanation  is 
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omitted  except  in  calling  attention  to  specific  considerations.     Tables 

21    and     22     show  the  successive  stages  obtained  through  the  steps  pre- 
viously discussed.     The  problem  consists  of  3  origins  and  5  destina- 
tions.      Cost  data  and  quantities  are  given  in  Table  21.     The  allocation 
to  routes  was  made  by  starting  in  the  northwest  corner  and  proceeding 
across  and  down.     The  requirements  are  fulfilled  by  six  cell  values, 
excluding  the  zero.     However,   in  this  method  the  number  of  vector  ele- 
ments in  the  basis  for  an    m  by  n    matrix  is   (m  •*-  n  -  l)    .     Therefore, 
7  vector  elements  must  be  used  (3+*  5  -  l)   •     The  optimal  program  may 
consist  of  less  than  (m  +- n  -  l)  routes  as  in  this  example  j   however,  dur- 
ing the  iteration  process  all  possible  basis  vectors  must  be  considered 
at  each  stage.     To  fulfill  the  requirement  of  7  vector  elements  in  the 
basis,   zero  is  arbitrarily  assigned  to     O2  Dp  •     The  zero  could  have  been 
assigned  to  any  of  the  9  remaining  cells  available,  and  this  would  re- 
quire different  row  and  column  numbers  than  given  in  Table  21.     The  trans- 
portation method  is  carried  out,     ?2\4     removed  and     P2£     brought  in.     The 
optimal  program  is  then  given  in  Table  22. 

There  are  5  non-stepping  stone  cells  with     z^-Cji     equal  to  zero, 
cells  13,  21,  31,   32,  and  33;  therefore,  those  cells  represent  alternate 
available  routes  at  the  same  cost  as  the  program  of  Table  22,  $230. 
V.'hen  there  is  more  than  one  basic  solution  the  other  solutions  consist 
of  a  convex  linear  combination  of  the  basic  solution.     Because  the  cost 
data  is  a  matrix,  if  a  number,  positive  or  negative,  is  added  to  any  en- 
tire row  or  to  any  entire  column,  the  new  problem  is  equivalent  to  the 
old  for  an  optimal  program."     By  subtracting  3>  h,  5»  6,  and  2  from  col- 
umns 1,   2,  3,  U,  £  respectively,  adding  2  to  row  1  and  adding  1  to  row 
2,  the  "cost"  data  of  Table  23  is  obtained.     The  route  can  be  programmed 
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according  to  the  new  relative  "costs";  the  -  5   "cost,"  cell  2>h,  may  be 
viewed  as  a  subsidy  route  with  the  maximum  amount  to  be  shipped  via  that 
route.     The  +  5>   "cost"  of  cell  15  is  avoided  by  use  of  the  many  alternate 
routes  of  0  "cost."     The  program  may  be  set  up  to  meet  all  requirements 
and  the  total  cost  obtained  through  the  use  of  the  original  cost  data. 

The  transportation  method  permits  cost  analysis  similar  to  the  sim- 
plex method,  for  factors  such  as:     what  cost  increase  is  permissible  for 
some  particular  route  to  meet  competitive  conditions  and  still  maintain 
an  optimum  program;  where  to  ship,  produce,  store;  and  similar  functions. 
The  objective  of  maximizing  some  functional  may  also  be  used  in  the 
transportation  method.     Symonds    '   discusses  an  oil  refinery  problem  of 
allocating  crude  oils  to  several  refineries  in  order  to  obtain  maximum 
possible  profit.     The  problem,  a  5  X  10  matrix,  is  solved  by  h  itera- 
tions. 

In  the  two  previous  illustrations,  a  separate  table  was  used  in 
each  iteration  for  clarity.     In  practice,  all  iterations  may  be  carried 
out  using  one  table  by  crossing  out  old  row  and  column  numbers  and  cell 
elements  changed.     It  is  a  characteristic  of  this  method  that  there  will 
be  no  undue  buildup  of  changed  values  in  any  one  cell.     Furthermore, 
there  is  no  need  of  rechecking  calculations  for  errors  as  a  positive 
check  is  provided  on  completion  of  each  iteration  by  the  use  of  new  row 
and  column  values. 
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ASSIGNMENT  METHOD 

Example  of  the  Assignment  Method 
The  assignment  method  is  a  special  case  of  the  transportation 
method  and  is  illustrated  in  a  problem  of  assigning  four  men  to  four 
jobs.  It  is  assumed  that  each  man  is  capable,  to  some  extent,  to  per- 
form each  job.  Furthermore,  it  is  assumed  that  the  degree  of  competence 
can  be  expressed  by  some  numeric  value,  such  as  a  score.  The  men  and 
corresponding  job  competence  in  each  job  is  indicated  in  Table  2U.  The 
objective  is  a  balanced  job  assignment;  it  is  desired  to  obtain  a  maxi- 
mum total  score  over  all  jobs,  where  a  high  score  indicates  competence. 
The  method  may  be  used  to  maximize  or  minimize  as  in  the  simplex  method. 
The  method  consists  in  "shipping"  fractions  of  men  to  the  various 

jobs;  mathematically  the  problem  is: 
m   n 


Maximize  (or  minimize)  *?     -*?   o^ix.^* 

i=l  j=l 
m 
Subject  to  *L   x.  •  =  total  man  required  on  the  j^*1  job  =  1  man 
i=l   J 

n 
£  x^  :  total  amount  of  i^*1    man  assigned  =  1  man 

Xj_^  z  amount  of  1th  man  assigned  to  the  j^*1  job. 
c..  z   score  of  ith  man  on  jth  job. 

Xy>0 

The  assignment  method  is  possible  because  of  the  following  theorem: 

Theorem:6     In  the  transportation  problem  when  marginal  data  are 
integers,  every  basic  solution  comes  out  in  integers. 

In  view  of  the  above,     Xjj     must  be  0  or  1  in  each  basic  solution  of  the 

problem.     To  minimize  total  score,  Table  2k  may  be  used  directly;  to 
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maximize  the  total  score,  the  scores  must  be  changed  to  negative  quan- 
tities. Thus  maximizing  the  original  score  is  equivalent  to  minimizing 
the  negative  of  the  score. 

To  maximize  total  score,  Table  2£  is  obtained  by  changing  the  scores 
of  Table  2k   to  negative  quantities.  Table  26  is  obtained  by  subtracting 
the  lowest  number  (-  97),  from  every  number  in  Table  25,  resulting  in 
all  positive  scores. 

For  each  iteration  process,  draw  the  minimum  number  of  horizontal 
or  vertical  lines  which  will  cross  out  all  zeros  in  the  array.  Look 
for  the  smallest  number  in  the  array  which  is  not  under  any  of  the  above 
lines.  Subtract  this  number  from  all  numbers  not  crossed  out  by  the 
horizontal  or  vertical  lines.  The  process  is  repeated  until  there  is  at 
least  one  zero  in  every  row  and  column.  This  condition  will  occur  when 
the  minimum  number  of  horizontal  and  vertical  cross  out  lines  is  equal 
to  the  number  of  rows  or  columns.  From  these  zeros,  select  the  zero  so 
that  it  is  the  only  one  in  each  row  and  column.  The  optimal  program  is 
given  in  Table  31  as  marked  by  the  asterisks:  man  1  to  job  h,   man  2  to 
job  3,  man  3  to  2,  and  man  h   to  job  1,  total  score  is  equal  to  339, 
(97  +  63  +•  89  +■  90). 

The  above  solution  requires  5  iterations  from  Table  26.  By  utiliz- 
ing the  matrix  characteristics  of  the  array  of  Table  26,  the  answer  may 
be  obtained  in  1  iteration.  Any  number  may  be  added  or  subtracted  to 
an  entire  row  or  an  entire  column  in  a  matrix.  Therefore,  Table  32  is 
obtained  directly  from  Table  26  by  subtracting  18  from  row  2,  3  from 
row  3,  and  7  from  row  h*     Table  33  is  obtained  by  subtracting  5  from 
column  2;   no  further  zeros  can  be  obtained  in  this  manner.  The  assign- 
ment method,  minimum  number  of  horizontal  or  vertical  lines  to  cross 
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out  all  zeros,  is  applied,  16  subtracted  from  the  remaining  array,  and 
the  optimal  program  of  Table  3h   obtained.  The  latter  is  identical  to 
Table  31. 

The  methods  discussed  work  only  for  a  square  matrix.  However,  the 
procedure  can  be  used  when  the  original  matrix  is  not  square  by  use  of 
fictitious  values  and  corresponding  high  or  low  "penalty"  values  which 
permit  solution  as  described  in  the  transportation  method.  In  addition, 
the  numerical  score  of  this  problem  could  have  been  replaced  by  time 
per  unit,  cost  per  unit,  quality  scale  per  vendor,  etc.,  depending  upon 
the  application. 
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Table  2h 

Jl     J2     J3  JU 

Mx    31    3U     73  97 

M2     55     50    63  79 

M3    37     89     9h  32 

M^     90    57     U2  37 


Table  25 

Jl      J2       J3  Jii 

31    -3U    -73  -97 

.55    -5o    -63  -79 

37    -89    -9U  -32 

-90    -57    -U2  -37 


Table  26 

Jl     J2  J3     JU 

-6^63  dlt— 9- 

U2     U7  3U     18 

60      8  3    65 

7    20  55    60 


Table  27 

Mx    66    63     2U 0- 

M2    39    UU    31    15 

M3    £? — 5 0    62 

M^      U    17    52    57 


Table  28 

66     63     2lt 0- 

35     U0     27     11 

-57 5 0-42- 

0    13    U8    53 


Table  29 
65    63 
2x    29 
Si      5 


3    6 


0    13    It  3    5 


Table  30 
Mx    66     58     2J| 
M2     2it     2U     16 
M3  -£? 0 9—* 

Mi,  — e — &— b8- 


Table  31 


50    lt2 
8       8 

-5f—o^ 


0*    8    U3    53 


B 

r 

3— 6J2 


Table  32 

66     63     2lt       0 

2U     29     36       0 

57      5       0    62 

0     13    U8     53 


Table 

33 

MX 

66     58 

2U 

0 

M2 

2it     2k 

16 

0 

M3 

n       P 

u    t 
lift     r 

UO     Jj 

r 

Table  3U 

5^    U2      |b 

8      8 
57 — 0*    b    6 


18    53 


DUAL  THEOREM 
The   "dual"  problem,  previously  mentioned,   is  developed  in  refer- 
ences  (l),   (U)  and  (6).     The  geometric  concept  of  the  direct  problem  and 
its  dual  problem  is  illustrated  for  a  two-dimensional  problem  in  ref- 
erence (21).     The  value  of  the  dual  problem  is  centered  on  the  theorem 
that  if  a  solution  exists  to  either  the  direct  problem  or  its  dual,  then 
a  solution  to  the  other  also  exists.     In  the  notation  previously  used, 
the  theorem  states: 

max.   £  CjXj  =  min.  2_  Wjbj_ 
3  i 

Thus,   in  formulating  the  model,   it  is  possible  to  work  with  either 
the  direct  problem  or  its  dual,  and  the  solution  of  either  automatically 
provides  the  solution  of  the  other.     The  answer  to  the  dual,  the    w^'s   , 
may  be  read  directly  from  the  final  tableau  from  the     Z-s-C-j     row  in  the 
columns  associated  with  the  initial  slack  vectors.     The     w^  ,  an  "oppor- 
tunity" value,   in  the  dual  problem  provides  information  on  the  rate  of 
change  of  the  functional  with  changes  to  the  capacities  or  constraints 
of  the  problem.     There  is  no  method  of  determining  wliich  form  of  the 
problem  will  yield  the  easiest  solution.     It  is  recommended  that  both 
forms  be  examined  before  attempting  to  solve  either  form. 

Use  of  the  dual  theorem  is  made  in  reference  (22),  a  warehousing 
problem  with  multiple  products  and  varying  prices.     It  also  illustrates 
special  computational  features  possible  due  to  the  structure  of  that 
problem. 
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SUMMARY 

Two  important  phases  of  linear  programming  are  formulation  of  the 
mathematical  model  and  the  computational  aspects. 

The  majority  of  linear  program  examples  published  to  date  involve 
one  or  more  of  the  methods  presented.  It  is  to  the  advantage  of  manage- 
ment and  industrial  engineers  to  become  acquainted  with  the  basic  methods 
and  capabilities  of  mathematical  programming,  particularly  in  industries 
that  do  not  maintain  a  staff  versed  in  this  field.  By  doing  so,  some  of 
the  necessary  groundwork  can  be  handled  by  present  personnel  to  reduce 
the  time  and  expense  in  utilizing  consultants  when  required.  Formulat- 
ing the  mathematical  model  requires  personnel  who  are  thoroughly  quali- 
fied in  the  use  of  the  various  techniques,  and  this  phase  determines  the 
success  or  failure  of  the  whole  process. 

For  a  large  problem,  solution  by  computer  is  the  only  efficient 
means,  not  because  of  the  complexity  of  the  mathematics  but  due  to  the 
quantity  involved.  A  transportation  problem  of  25  origins  and  60  des- 
tinations required  nine  man  days  by  hand  computation  techniques.  An 
IBM  card  programmed  calculator  can  solve  a  20  by  20  matrix  in  about 
three  hours.  The  Rand  Corporation  (23)  simplex  code  can  solve  100 
equations  in  about  five  hours,  and  their  code  for  a  transportation  type 
problem  can  solve  a  300  equation  system  in  about  one  hour. 

Rental  time  at  data  processing  centers,  located  in  various  parts  of 
the  country,  offer  one  possibility  for  snail  companies  unable  to  justify 
the  large  capital  Investment  of  a  computer. 

References  (2U)  and  (25)  discuss  methods  of  machine  matrix  inver- 
sion. The  method  of  reference  (2U)  is  not  well  suited  to  desk  calculator 
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application  because  of  the  amount  of  record  keeping  at  each  step.  Ref- 
erence (26)  provides  a  complete  description  on  setting  up  a  problem  for 
solution  by  an  IBM  701  computer. 
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